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24 Abstract 

25 The development of porcine epidemic diarrhea virus (PEDV) antibody-based assays is important 

26 for detecting infected animals, confirming previous virus exposure, and monitoring sow herd 

27 immunity. However, the potential cross-reactivity among porcine coronaviruses is a major 

28 concern for the development of pathogen-specific assays. In this study, we used serum samples 

29 (n = 792) from pigs of precisely known infection status and a multiplex fluorescent microbead- 

30 based immunoassay and/or enzyme-linked immunoassay platform to characterize the antibody 

31 response against PEDV whole-virus (WV) particles and recombinant polypeptides derived from 

32 the four PEDV structural proteins, i.e., spike (S), nucleocapsid (N), membrane (M), and envelope 

33 (E). Antibody assay cut-off values were selected to provide 100% diagnostic specificity for each 

34 target. The earliest IgG antibody response was observed at days 7-10 post-infection, mainly 

35 directed against SI polypeptides. With the exception of non-reactive protein E, we observed a 

36 similar antibody ontogeny and pattern of seroconversion for SI, N, M, and WV antigens. 

37 Recombinant SI provided the best diagnostic sensitivity, regardless of PEDV strain, with no 

38 cross-reactivity detected against transmissible gastroenteritis virus (TGEV), porcine respiratory 

39 coronavirus (PRCV), or porcine deltacoronavirus (PDCoV) pig antisera. The WV particles 

40 showed some cross-reactivity against TGEV Miller and TGEV Purdue antisera, while N protein 

41 presented some cross-reactivity against TGEV Miller. The M protein was highly cross-reactive 

42 against TGEV and PRCV antisera. Differences in the antibody response against specific PEDV 

43 structural proteins have important implications in the development and performance of antibody 

44 assays for the diagnosis of PEDV enteric disease. 

45 Key words: PEDV; recombinant structural proteins; whole virus; multiplex FMIA; ELISA, 
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47 1. Introduction 

48 Porcine epidemic diarrhea vims (PEDV) is an enveloped, singled-stranded, positive- 

49 sense RNA vims that belongs to the order Nidovirales, family Coronaviridae, subfamily 

50 Coronavirinae, and genus Alphacoronavirus (1). The PEDV genome (~28 kb) consists of seven 

51 open reading frames (ORFs) (2). The 5’ two-thirds contain the replicase-transcriptase ORF1 

52 (overlapping ORFla and ORFlb), followed by five ORFs encoding four structural proteins and 

53 one strain-specific accessory protein in the following order: spike (S), ORF3 (accessory), 

54 envelope (E), membrane (M), and nucleocapsid (N) (3). 

55 PEDV was first reported in Europe in the early 1970s as the causative agent of PED (7). 

56 PEDV classical CV777-like strains were subsequently reported in Europe and Asia (8), but 

57 PEDV was absent from the Americas, Africa, and Oceania prior to 2013 (9). The emergence of 

58 high-virulence PEDV strains was first recognized in late 2010 in China, with outbreaks reported 

59 in April 2013 in U.S. (5). Subsequently, high-virulence PEDV strains have been the cause of 

60 major economic loss in the swine industry worldwide, producing high mortality in neonatal 

61 piglets and high morbidity, but moderate mortality, in weaned pigs (4-6). The emergent PEDV 

62 strains are genetically distinct from the classical CV777-likes strains that continue to circulate in 

63 the field (5, 10, 11). Based on differences in the S gene and virulence, emerging PEDV strains 

64 can be divided into non S-INDEL (S gene insertions and deletions) and S-INDEL strains (9, 12). 

65 Overall, S-INDEL strains cause lower mortality than the high-virulence non-S INDEL strains 

66 (13, 14). 

67 In addition to PEDV, three other porcine enteric coronavimses have been described: 

68 transmissible gastroenteritis coronavims (TGEV) (15), porcine deltacoronavirus (PDCoV) (16), 

69 and a recently described swine enteric coronavirus (SeCoV) that emerged by recombination 
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70 between TGEV and PEDV (17). Enteric coronaviruses primarily infect villous enterocytes, 

71 causing atrophic enteritis that leads to malabsorptive diarrhea (4, 5, 18). In general, PEDV and 

72 TGEV are considered more virulent than PDCoV, but the three pathogens are clinically and 

73 histopathologically indistinguishable (5, 14, 19). Porcine respiratory coronavirus (PRCV) has a 

74 predilection for the respiratory tract, but PRCV is an S gene deletion mutant of TGEV and 

75 remains on the list of enteric coronavirus differentials. 

76 The differential diagnosis of porcine enteric coronaviruses relies on laboratory direct- 

77 detection methods, e.g., polymerase chain reaction (PCR) methods, immunohistochemistry, 

78 fluorescent in situ hybridization, and direct immunofluorescence in tissues (20-23). Antibody- 

79 based assays play an important role in detecting infection and evaluating immunity, but antibody 

80 cross-reactivity between porcine enteric coronaviruses is a major concern. As part of the process 

81 of developing PEDV-specific antibody assays, we experimentally inoculated pigs with each of 

82 the porcine coronaviruses (PEDV, TGEV, PRCV, and PDCoV) and characterized the antibody 

83 response to recombinant polypeptides derived from PEDV structural proteins, and to the intact 

84 PEDV virion using a multiplex fluorescent microbead-based immunoassay (FMIA) and a whole- 

85 virus (WV) ELISA. The final aim of this project was to identify highly sensitive and specific 

86 PEDV antigen targets for the antibody-based differential diagnosis of coronavirus-related enteric 

87 disease. 

88 

89 2. Material and Methods 

90 2.1. Experimental design 

91 The study was conducted under the approval of the Iowa State University Office for 

92 Responsible Research. Seven-week-old pigs (n = 72) were purchased from a conventional wean- 
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93 to-finish farm with no previous history of porcine coronavirus infections. The pigs were 

94 prescreened for evidence of infection with PEDV, TGEV, PRCV, and PDCoV. Pig fecal swabs 

95 were tested by a PEDV N gene-based rRT-PCR (21) and PDCoV M gene-based rRT-PCR (14), 

96 while pig fecal and nasal swabs were tested by TGEV (S gene)/PRCV (N gene)-based 

97 differential rRT-PCR (20). The pigs’ serum samples were tested with the PEDV 

98 immunofluorescence assay (IFA) (21), PEDV WV ELISA (24), TGEV/PRCV differential ELISA 

99 (Svanova, Sweden), and PDCoV IFA (14). Animals (n = 72) were randomized into six groups; 

100 each group consisted of 12 pigs in one room, with 6 pens per room and 2 pigs per pen. Details 

101 related to virus strains and the routes of experimental inoculation are presented in Table 1. The 

102 age-related infectious dose for each virus was previously determined in a pilot study (data not 

103 shown). The pigs were closely observed twice daily for clinical signs throughout the study. A 

104 total of 792 serum samples were collected from each group on day post-infection (DPI) -7, 0, 3, 

105 7, 10, 14, 17, 21, 28, 35, and 42. Virus shedding within groups and absence of cross- 

106 contamination between groups during the observation period (-7 to 42 DPI) was confirmed by 

107 rRT-PCR, whereas seroconversion within inoculation groups was confirmed by ELISA or IFA 

108 (data not shown). At 42 DPI, all pigs were humanely euthanized by penetrating captive bolt 

109 (Accles and Shelvoke, Ltd., Sutton Coldfield, UK) followed by exsanguination. 

110 Five different polypeptides corresponding to the four PEDV structural proteins (S, N, M, and 

111 E) were recombinantly generated; PEDV WV particles were purified from cell culture by 

112 ultracentrifugation, as described elsewhere (24). The PEDV recombinant and WV target antigens 

113 were used to develop a multiplex (6-plex) FMIA platform to analyze the antibody response to 

114 each target using sera from pigs experimentally inoculated with PEDV, TGEV Miller and Purdue 

115 strains, PRCV, PDCoV, or sham control. The diagnostic sensitivity and specificity and the 
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116 analytical specificity were evaluated for each individual antigen. 

117 2.2. Generation of PEDV recombinant spike (rS)-derived proteins 

118 The coding region of the amino-terminal receptor-binding (SI) domain derived from a 

119 consensus sequence (Figure 1) based on ten PEDV non S-INDEL strains (2,151 nt) or five S- 

120 INDEL strains (2,142 nt) (Table 2), respectively, were codon-optimized for expression in 

121 mammalian cells and synthetically produced (Shanghai Genery Biotech Co., Ltd.) with the 

122 addition of a 5’ terminal eukaryotic native signal (encoding MKSLTYFWLFLP VLSTLSLP) and 

123 a 3’ terminal Tobacco Etch Vims (TEV) cysteine protease site (encoding ENLYFQS), followed 

124 by the Fc portion of human IgGl (Genbank #JX292764.2). The genes were amplified or 

125 manipulated using the forward primers F1 -(5 ’-TAA ACG GAT CTC TAG CGA ATT CGC CGC 

126 CAC CAT GAA GAG CCTG-3’) and F2-(5’-CTT CCA GAG CGG CTC CGA CAA GAC CCA 

127 CAC CGT CGA GTG CCC ACC GTG CCC AG-3’), and reverse primers Rl-(5’-CGA GCG 

128 GCC GCT AGA AGC TTT CAT TTA CCC GGA GAC AGG GAG-3’) and R2-(5’-GTC GGA 

129 GCC GCT CTG GAA GTA CAG GTT CTC GTG ATA GAA GAA TCC GGG CAG-3’). 

130 Amplicons were cloned into pNPM5 eukaryotic expression vector (Novoprotein, Short Hills, NJ, 

131 USA) and the recombinant plasmids were transiently transfected into human embryonic kidney 

132 (F1EK) 293 cells (lxl O 6 cells/ml) (Invitrogen, Thermo Fisher Scientific, Grand Island, NY, USA) 

133 using polyethylenimine (PEI) (Thermo Fisher Scientific) at an optimal 1:4 ratio (plasmid:PEI, 

134 w/w). Transfected F1EK293 cells were grown in serum-free FreeStyle™ 293 Expression Medium 

135 (Gibco®, Life Technologies, Carlsbad, CA, USA) at 37°C with 5% CO 2 by orbital shaking at 

136 120 rpm. At day 5 post-transfection, culture supernatants were harvested by centrifugation at 

137 3,500 xg for 20 min, and filter-sterilized (0.45 pm filter). The soluble expression of Fc-tag fused 

138 SI proteins (107.5 kDa) was confirmed by 12% dodecyl sulfate-polyacrylamide gel 
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139 electrophoresis (SDS-PAGE). PEDV Sl-Fc proteins were enzymatically cleaved by incubation 

140 with TEV (20 IU/mg sample) for 3 h at 25°C under endotoxin control, and purified from culture 

141 supernatant by protein A chromatography (GE Healthcare, Pittsburgh, PA, USA), and nickel 

142 (Ni)-chelating Sepharose Fast Flow (SFF) affinity chromatography (GE Healthcare), according 

143 to the manufacturer’s instructions. Purified l'Sl non S-INDEL (717 aa) and l'Sl S-INDEL (714 

144 aa) proteins were dialyzed against phosphate-buffered saline (PBS) (10 mM phosphate and 150 

145 mM NaCl, pH 7.4) and analyzed by SDS-PAGE and Western blot. 

146 2.3. Generation of PEDV recombinant nucleocapsid (rN) protein 

147 An Escherichia coli ( E. coli)- codon optimized consensus version (Figure 1) obtained from a 

148 multiple sequence alignment (non S-INDEL, S-INDEL, and CV777 strains) of the full-length 

149 PEDV N (1,356 nt) gene (Table 2) was synthesized in vitro (Shanghai Genery Biotech Co., Ltd., 

150 Shanghai, China). The gene was amplified using the forward primer (5’-CAT CAT CAT CAT 

151 CAT CAT ATG GCA TCT GTT AGC TTT CAG GAT CG-3’) and reverse primer (5’-AGA CTG 

152 CAG GTC GAC AAG CTT TTA ATT GCC GGT ATC GAA GAT C-3’). The amplicon was 

153 cloned into pCold II expression plasmid (Novoprotein Scientific Inc., Shanghai, China), 

154 confirmed by sequencing (Genewiz Inc., Suzhou, China), and then transformed into the E. coli 

155 BL21 (DE3) host strain (Invitrogen™, Carlsbad, CA, USA). The transformants were 

156 resuspended and grown in 1 liter of Luria-Bertani (LB) medium (Invitrogen™) containing 100 

157 ug/ml ampicillin, at 16°C by shaking at 250 rpm. When an Agoo of 0.9 was reached, 0.1 mM 

158 isopropyl—thio-D galactopyranoside (IPTG) was added, and cultures were grown for an 

159 additional 16 h at 16°C. Cells were chilled at 4°C and harvested by centrifugation at 3,500 *g for 

160 15 min, resuspended in 150 ml of lysis buffer (20 mM Tris-HCl and 500 mM NaCl, pH 8.0), and 

161 disrupted by ultrasonication (Vibra-Cell 1M sonicator, Sonics & Materials, New Town, CT, USA). 
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162 The crade extracts were centrifuged at 15,000 xg for 30 min at 4°C and the soluble expression of 

163 histidine (his)-tag fused N-PEDV protein (50.4 kDa) was confirmed by SDS-PAGE analysis. The 

164 rN protein (452 aa) was purified from the soluble fraction by Ni-chelating SFF affinity 

165 chromatography (GE Flealthcare) according to the manufacturer’s instructions. After separation 

166 of thrombin cleavage products on the Nr column, FliTrap™ Phenyl High Performance (TIP) 

167 (GE Flealthcare) hydrophobic interaction chromatography followed by a FliTrap™ SulfoPropyl 

168 (SP) TIP strong cation exchange chromatography (GE Flealthcare) were consecutively applied 

169 according to the manufacturer’s instructions. Protein elutions were pooled, dialyzed against PBS, 

170 pFI 7.4, at 4 °C, and analyzed by SDS-PAGE and Western blot. 

171 2.4. Generation of PEDV recombinant membrane (rM) protein 

172 An E. coli- codon optimized consensus version (Figure 1) from multiple alignment (non S- 

173 INDEL, S-1NDEL, and CV777 strains) of the PEDV M (429 nt) gene encoding a truncated 

174 fragment corresponding to the C-terminal intra-virion topological domain of the M protein 

175 (Table 2) was synthetically generated (Shanghai Genery Biotech Co., Ltd.). The gene was 

176 amplified using the forward primer (5’-CAT CAT CAT CAT CAT CAT ATG TTT GTG AAT 

177 AGT ATT CGC TTATGG-3’) and reverse primer (5’-AGACTG CAG GTC GAC AAG CTT 

178 TTA A AC CAG ATG CAG AAC TTT TTC G-3’). PCR products were cloned, sequenced, and 

179 transformed into BL21 (DE3) cells as described above for the N protein. The truncated his- 

180 tagged fusion M-PEDV polypeptide (15.7 kDa; 143 aa), expressed in the precipitate of cell lysate 

181 as an inclusion body, was solubilized using a denaturing buffer (20 mM Tris-FICl, 300 mM NaCl, 

182 8 M urea, and 5 mM P-mercaptoethanol; pFI 8.0), and purified by Ni-chelating chromatography 

183 (GE Flealthcare). Elutions were pooled, dialyzed against 20 mM phosphate buffer, 150 mM 

184 NaCl, and 6 M urea, pFI 7.4, at 4 °C, and analyzed by SDS-PAGE and Western blot. 
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185 2.5. Generation of recombinant envelope (rE) small membrane protein 

186 A commercial E. co/Z-expressed PEDV (strain CV777) recombinant his-tagged E protein (76 

187 aa) was purchased from Cusabio (CSB-EP771125PPW; Cusabio Biotech Co., MD, USA) and it 

188 sequence is presented in Figure 1. 

189 2.6. PEDV multiplex fluorescent microbead-based immunoassay (FMIA) 

190 The covalent coupling of purified PEDV WV antigen and recombinant polypeptides to 

191 microbead sets was performed as previously described (25). A total of 25 pg of each protein, i.e., 

192 rSl non S-INDEL (0.6 mg/ml; bead region 53), rSl S-INDEL (0.3 mg/ml; bead region 45), rN 

193 (0.44 mg/ml; bead region 64), rM (1 mg/ml; bead region 12), and rE (1 mg/ml; bead region 54), 

194 and 60 pi of purified WV antigen were coupled to 5x 10 6 carboxylated magnetic microspheres 

195 (MagPlex®-C Microspheres, Luminex Corp., Austin, TX, USA). Each set of coupled beads were 

196 individually evaluated using different assay buffers to exclude matrix inhibitory effects, and to 

197 select optimum buffer combination working for all bead sets. After verified that each bead sets 

198 worked properly in a single-analyte assay, they were sequentially added into progressively larger 

199 multiplex assay to exclude interference between bead regions. Serum samples were diluted to 

200 1/50 in assay buffer (Sea Block Blocking Buffer, Thermo Fisher Scientific) and mixed with 50 pi 

201 of the bead suspension (~2,500 beads per well for each target) in each well (Bio-Plex Pro™ flat 

202 bottom plates, Bio-Rad Laboratories Inc., Flercules, CA, USA). The plates were incubated for 30 

203 min on a microplate shaker (VWR®, Radnor, PA, USA) at 400 rpm, and washed three times 

204 with 200 pi of 0.1 M PBS, pFI 7.4, containing 0.1% Tween-20 (PBST). All incubations were 

205 performed at ~22°C in a dark environment. Next, 50 pi of biotin-labelled goat anti-pig IgG (Fc) 

206 (Bethyl Laboratories Inc., Montgomery, TX, USA) at 1/3,000 in assay buffer was added to each 

207 well, followed by a 30 min incubation. After a washing step, 50 pi streptavidin phycoerythrin 

208 (SAPE; Moss Inc., Pasadena, MD, USA) at 1/100 in assay buffer was added to each well and the 
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209 plates were incubated for 30 min. After an additional wash step, the microspheres were 

210 resuspended in 100 pi of assay buffer and analyzed using a dual-laser Bio-Plex® 200 instrument 

211 (Bio-Rad Laboratories, Inc.). The events were gated to exclude doubles and other aggregates 

212 (Bio-Plex Manager™ software 6.0, Bio-Rad Laboratories, Inc.). Antibody response, reported as 

213 median fluorescent intensity (MFI), was expressed as the sample-to-positive (S/P) ratio: 

S/P ratio - ( sam Pl e FI ~ background control mean FI) _ 

(positive control mean FI - background mean FI) 

214 

215 2.7. Data analysis 

216 The cutoff and diagnostic performance of each individual PEDV antigen was determined by 

217 receiver operating characteristic (ROC) analysis (SAS® Version 9.4, SAS® Institute, Inc., Cary, 

218 NC, USA) using 6-plex FMIA and/or WV ELISA test results. Porcine coronavirus-negative 

219 samples (n = 252) were used to estimate diagnostic specificity and PEDV-positive samples (n = 

220 72) collected from the 12 PEDV-inoculated pigs between DPI 14 to 42 were used to estimate 

221 diagnostic sensitivity, time of detection, and over-time detection through the observational 

222 period. The analytical specificity (cross-reactivity) of each PEDV antigen was evaluated using 

223 samples (n = 384) collected between DPI 7 and 42 from animals inoculated with TGEV Miller 

224 and Purdue strains, PRCV, and PDCoV. 

225 Pearson’s chi-square test was used to detect differences in the number of seropositive 

226 animals over time (by DPI) among the individual antigens. Fisher’s exact test was used when 

227 25% of the cells in a contingency table contained counts of <5 (SAS® Version 9.4, SAS® 

228 Institute, Inc., Cary, NC, USA). One-way ANOVA with Tukey’s correction was used for multiple 

229 comparisons with alpha=0.05 (GraphPad Prism® 6). Specifically, we compare the antibody 

230 response between inoculation groups by day post inoculation for each antigen target. 
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231 3. Results 

232 3.1. Dynamics of antibody responses to different PEDV antigens after experimental 

233 inoculation 

234 The IgG serum antibody responses to individual PEDV antigens (rSls [S-INDEL and non S- 

235 INDEL], rN, rM, rE, and WV) were evaluated over time (DPI -7 to 42) in pigs inoculated with 

236 PEDV, TGEV Miller, TGEV Purdue, PRCV, PDCoV, or negative control by 6-plex FMIA 

237 (Figure 2) or PEDV WV ELISA (Figure 4A). In the PEDV-inoculated group, similar antibody 

238 dynamics against rSls, rN, rM, and WV antigens were observed, with significantly higher 

239 (p<0.05) antibody levels than the negative control group at DPI >10. Likewise, antibody 

240 responses against rSls (FMIA), rN (FMIA), and WV (ELISA) in the PEDV-inoculated group 

241 were significantly higher (p < 0.05) at DPI >10 than the TGEV, PRCV, and PDCoV inoculation 

242 groups. For rM protein, there were significant differences (p < 0.05) between the PEDV and 

243 PDCoV inoculation groups; however, no significant differences in rM response were found in 

244 pairwise comparisons between the following treatment groups: PEDV vs TGEV Miller (over 

245 time), PEDV vs TGEV Purdue (DPI -7 to 10 and DPI 28 to 42), and PEDV vs PRCV (DPI -7 to 

246 10 and DPI 42). No differences in the FMIA IgG S/P ratio antibody response to SI from the S- 

247 INDEL strain compared to SI from the non S-INDEL strain were observed over time in any 

248 inoculation group. 

249 3.2. Diagnostic sensitivity and specificity of PEDV antigens 

250 The optimal S/P cut-off points were selected to ensure 100% diagnostic specificity. Selected 

251 cut-off values and diagnostic sensitivity associated for each target tested by FMIA (rS 1 non S- 

252 INDEL, rS 1 S-INDEL, rN, rM, rE, and WV) or ELISA (WV) is presented in Table 1. 

253 The time of detection and the over-time detection of PEDV antibodies to different PEDV 

254 antigens is presented in Table 2. The first PEDV antibody-positive pigs were observed on DPI 7 
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255 for all PEDV antigens, except rE protein. Recombinant E protein showed poor antigenicity, 

256 inducing seroconversion in only one animal infected with PEDV (1/12) at DPI 35. With the 

257 exception of the non-reactive rE protein, no significant differences (p > 0.05) in the number of 

258 seropositive animals were observed over time among the PEDV antigens. Likewise, no 

259 difference was observed over time in the proportion of seropositive animals detected by both 

260 rSls derived from the PEDV non S-INDEL and S-INDEL strains. Interestingly, the rate of 

261 seropositivity using the WV antigen was higher when it was used on the ELISA platform as 

262 opposed to the FMIA platform. 

263 3.3. Cross-reactivity of PEDV antigens against PEDV-related porcine enteric coronaviruses 

264 The cumulative distribution of serum antibody responses (S/P) against specific PEDV 

265 markers is given in Figure 3 (6-plex FMIA) and Figure 4B (WV ELISA) for each inoculation 

266 group. Based on the selected cut-off points, the overall analytical specificity of individual PEDV 

267 antigens was 100% for both rSl polypeptides and rE protein, 98.7% for rN protein, 51.6% for rM 

268 polypeptide, 95.6% for WV in the 6-plex FMIA, and 96.3% for WV ELISA. The cross-reactivity 

269 of individual antigens against antiserum obtained over time from pigs inoculated with PEDV- 

270 related porcine coronaviruses (TGEV Miller, TGEV Purdue, PRCV, and PDCoV) tested by the 6- 

271 plex FMIA and/or WV ELISA is presented in Table 2. PEDV rN protein cross-reacted against 

272 TGEV Miller antisera (2/12 pigs) at DPI 7 and 14. Similarly, WV FMIA detected one positive 

273 pig from both the TGEV Miller and PRCV infected groups, and WV ELISA detected one 

274 positive animal from the TGEV Purdue infected group. The WV antibody response against 

275 PEDV was higher and discriminated better between groups with ELISA (Figure 4A) compared to 

276 FMIA (Figure 2). The rM protein was highly reactive against pigs experimentally inoculated 

277 with TGEV Miller (12/12), TGEV Purdue (9/12), and PRCV (7/12) between DPI 7-42. 
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278 Moreover, rM protein was the only antigen showing some cross-reactivity in pigs (1/12 pigs at 

279 DPI 35) inoculated with PDCoV. 

280 The comparison of the amino acid sequence homology of the 5 purified PEDV polypeptides 

281 recombinantly generated (i.e., SI non S-INDEL and S-INDEL, N, M, and E) with the 

282 homologous regions corresponded to porcine coronaviruses PEDV, TGEV Miller, TGEV Purdue, 

283 PRCV, and PDCoV used during experimental inoculation is presented in Table 3. 

284 

285 4. Discussion 

286 The pig immune system has the ability to recognize specific PEDV proteins and to 

287 respond by producing specific antibodies (21, 24, 26-28) and the development of PEDV antibody 

288 assays is important for detecting infection, confirming previous virus exposure, and monitoring 

289 levels of immunity. To characterize the antibody response against primary PEDV infection, five 

290 recombinant polypeptides derived from the four PEDV structural proteins (S, M, N, and E) and 

291 a WV antigen were generated. Proteins S and M were truncated based on gene size, 

292 hydrophobicity, and antigenic features to facilitate overexpression and efficient purification. A C- 

293 terminal intravirion topological domain of the M protein, previously identified as the major 

294 immunodominant region in the M protein, was selected (29). PEDV S protein can be structurally 

295 or functionally divided into two subunits: SI (N-terminal globular head), which binds to the host 

296 cell receptors, and S2 (C-terminal membrane-bound stalk), which is responsible for membrane 

297 fusion (30, 31). The SI subunit was selected because it contains major antigenic sites and 

298 antiviral neutralizing determinants in many animal coronaviruses, including TGEV (32, 33). 

299 Elowever, S1 diverges in sequence even among species of a single coronavirus, thereby 

300 contributing to the broad host range of coronaviruses, whereas the S2 subunit is the most 
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301 conserved region of the protein (34, 35). Therefore, two recombinant SI proteins derived from 

302 PEDV non S-INDEL and S-INDEL strains, respectively, were generated and comparatively 

303 evaluated. 

304 The analysis showed that serum antibody from PEDV-inoculated pigs reacted to all 

305 proteins evaluated, except rE. These data are in agreement with previous results on the utility of 

306 N (26), S (28), M (27), and WV antigen (24) for PEDV antibody detection. In general, the 

307 magnitude of an antibody response is dependent upon the amount of antigen presented to the 

308 immune system. The M protein is the most abundant protein in the virion envelope (36), the N 

309 protein is the most abundant coronavirus antigen produced throughout infection (37), and the S 

310 protein forms the crown-like projection of the viral surface (38). In contrast, E protein is present 

311 only in small amounts in infected cells and the viral envelope (39); which may explain its poor 

312 antigenicity. 

313 A difference in the WV antibody response was noted between the FMIA and ELISA 

314 formats. The lower response in the FMIA assay can be attributed to the fact that the Luminex® 

315 carboxylated magnetic microbeads bind WV (intact) PEDV particles less efficiently than they 

316 bind proteins or polypeptides with free amino groups. Likewise, the highly charged polystyrene 

317 microwells utilized on the ELISA platform bind WV PEDV particles most efficiently, providing 

318 better detection. 

319 PEDV IgG antibodies against SI (S-INDEL and non S-INDEL), N, M, and WV were 

320 first detected between 7 to 10 DPI and thereafter for the duration of the study. However, our data 

321 confirmed the antibody response to SI proteins as the most sensitive marker of early PEDV 

322 infection and the best choice for reliable detection of weak seroresponders (Table 2). Notably, 

323 antibodies produced against the PEDV non S-INDEL (USA/IN/2013/19338E) strain used for 
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324 experimental inoculation reacted similarly against the two SI proteins generated from a 

325 consensus between different S-INDEL and non S-INDEL PEDV strains, respectively. Previously, 

326 Chen et al. (2016) demonstrated that the antibodies against U.S. PEDV S-INDEL and non S- 

327 INDEL strains cross-reacted and cross-neutralized both strains in vitro (40). Our results 

328 confirmed that the high percentage of homology in the amino acid sequence of SI protein among 

329 PEDV non-S INDEL and S-INDEL strains (>90%) assures that SI could be used as a diagnostic 

330 marker for a wide range of PEDV strains. 

331 In the field, pigs are exposed to different coronaviruses that are known to share genetic 

332 and antigenic traits that may contribute to false-positive results. Recent evidence suggesting 

333 antibody cross-reactivity between PEDV and TGEV (13) and between PEDV and PDCoV (41) 

334 raised concerns about the specificity of PEDV serologic testing using assays based on targets 

335 containing the most conserved regions (e.g., N, M, and WV). Therefore, the serologic cross- 

336 reactivity between individual PEDV antigens and other swine coronaviruses was evaluated by 

337 testing antisera from pigs experimentally inoculated with TGEV (Miller and Purdue strains), 

338 PDCoV, and PRCV using 6-plex FMIA and WV ELISA. 

339 The SI-derived antigens proved to be the best candidate for antibody-based differential 

340 testing of porcine coronaviruses based on a complete absence of detectable cross reactivity 

341 (100% analytical specificity). Protein N showed some cross-reactivity (2/12 pigs) against TGEV 

342 Miller, but not against TGEV Purdue antisera, and only at the early stages post-infection (DPI 7 

343 to 14). This “transient” cross-reactivity may explain why it has been reported infrequently. In 

344 contrast, PEDV M protein showed the highest percentage of homology with other porcine 

345 coronaviruses (Table 4) and should be ruled out as a marker for differential diagnosis, as it was 

346 highly cross-reactive against TGEV strains and PRCV. 
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347 Differences in cross-reactivity among TGEV strains cannot be explained simply on the 

348 basis of amino acid sequence homology (Table 3). Lin et al. (2015) reported one-way cross- 

349 reactivity between different PEDV strains and TGEV Miller antisera, but not TGEV Purdue 

350 antisera (42). Further analysis identified at least one epitope on the N-terminal region of the 

351 PEDV/TGEV N protein that contributed to this cross-reactivity (13). In the present study, a full- 

352 length N protein was produced. Truncation of the N-terminal region could help avoid possible 

353 cross-reactivity (42, 43). Alternatively, cross-reactivity could be corrected by increasing the cut- 

354 off from S/P >0.59 to S/P >0.80, but this would reduce the diagnostic sensitivity from to 83.3% 

355 to 76.3%. Similarly, the cross-reactivity of the WV antigen against TGEV Miller (1/12 pigs) and 

356 TGEV Purdue (1/12) antisera revealed by ELISA (Table 2) could be corrected using a more 

357 conservative cut-off (S/P >1.3 instead of S/P >0.7) without severely impacting diagnostic 

358 sensitivity (from 88.8% to 83.3%). 

359 This study demonstrated that variations in the antibody response against different PEDV 

360 structural proteins may have important implications in the diagnosis of PEDV enteric disease. 

361 We also successfully identified targets of interest (e.g., SI) for the diagnosis of PEDV, providing 

362 a truly molecular immunological view of antigenic distribution and a complete antibody cross- 

363 reactivity profile between PEDV and other porcine enteric coronaviruses. 
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512 Tables 


513 


514 Table 1. Optimal S/P cut-off points and diagnostic sensitivity of serum IgG antibody 

515 response to each PEDV antigen by 6-plex FMIA or WV ELISA. 

516 


517 


PEDV marker 

Assay 

Cut-off No samples tested positive/Total No positive samples (%); 95% Cl 
(S/P)* 

SI (non S-INDEL) 


0.10 

69/72 (95.8); 88.4-98.6 

SI (S-INDEL) 


0.11 

68/72 (94.4); 86.5-97.8 

N 

6-plex 

0.59 

59/72 (81.9); 71.5-89.1 

M 

FMIA 

0.21 

58/72 (80.6); 69.9-88.0 

E 


0.90 

1/72(1.4); 0.25-7.5 

WV particles 


0.67 

48/72 (66.7); 55.2-76.5 

WV particles 

ELISA 

0.70 

63/72 (87.5); 77.9-93.3 


*The optimal S/P cut-off points were selected to ensure 100% diagnostic specificity 


518 


519 Table 2. Detection of serum IgG antibody response among inoculation groups by day post 

520 inoculation. Data are presented as number of positive samples (above selected S/P cut-off 

521 values) to each PEDV antigen by 6-plex FMIA or WV ELISA. Dark grey shaded boxes 

522 indicate true positive results and light grey shaded boxes indicate false positive results 

523 (cross-reactivity). 


Inoculation 

group 

No pigs 

PEDV marker 

Assay 

Number pigs positive by day post-inoculation 


-7 

0 

3 

7 

10 

14 

17 

21 

28 

35 

42 



SI (non S-INDEL) 


0 

0 

0 

6 

11 

11 

11 

12 

12 

12 

11 



SI (S-INDEL) 


0 

0 

0 

5 

10 

11 

11 

12 

12 

11 

11 



N 


0 

0 

0 

2 

6 

10 

10 

10 

10 

10 

9 

PEDV 

12 

M 

6-plex FMIA 

0 

0 

0 

1 

8 

10 

10 

10 

10 

9 

9 



E 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 



WV particles 


0 

0 

0 

0 

3 

7 

9 

9 

8 

8 

7 



WV particles 

Single ELISA 

0 

0 

0 

2 

7 

11 

10 

11 

11 

10 

10 



SI (non S-INDEL) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



SI (S-INDEL) 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



N 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

TGEV Purdue 

12 

M 

6-plex FMIA 

0 

0 

0 

3 

9 

9 

9 

9 

8 

8 

8 



E 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



WV particles 


0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 



WV particles 

Single ELISA 

0 

0 

0 

1 
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524 


525 Table 3. Amino acid sequence homology of PEDV recombinant polypeptides Spike (SI 

526 globular head), Nucleocapsid (N; full length), Membrane (M; C-terminal intravirion 

527 topological domain), and Envelope (E; full length) compared to homologous regions of 

528 PEDV-related porcine coronaviruses TGEV Miller (ATCC VR-1740), TGEV Purdue 

529 (ATCC VR-763), PRCV (ATCC VR-2384), PDCoV (USA/IL/2014) used during 

530 experimental inoculation. Sequence analysis was performed using Unipro UGene 

531 multiplatform software (Version 1.25, Novosibirsk, Russia). 

532 



% homology amino acid sequence of PEDV recombinant polypeptides 

Strains used for 
experimental inoculation 

SI PEDV 
non S-INDEL 
(consensus) 

SI PEDV 
S-INDEL 
(consensus) 

N PEDV 
(consensus) 

M PEDV 
(consensus) 

E PEDV 
(CV777 strain) 

PEDV non S-INDEL 

100% 

92% 

100% 

100% 

98% 

(USA/IN/2013/19338E) 

TGEV Purdue 

28% 

28% 

31% 

55% 

28% 

(ATCC VR-763) 

TGEV Miller 

28% 

28% 

31% 

55% 

28% 

(ATCC VR-1740) 

PRCV 

29% 

29% 

31% 

55% 

30% 


(ATCC VR-2384) 
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PDCoV 23% 23% 20% 20% 
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536 Figure legends 

537 Figure 1. Amino acid sequences of 5 recombinant polypeptides derived from the 4 PEDV 

538 structural proteins named: 1) Envelope (E) full length protein (CV777 PEDV strain); 2) Full 

539 length nucleocapsid (N) protein from a consensus among PEDV strains (CV777, non S-INDEL, 

540 and S-INDEL); 3) Truncated (C-terminal intravirion topological domain) membrane (M) protein 

541 from a consensus among PEDV strains (CV777, S-INDEL, non S-INDEL); 4) Truncated spike 

542 protein (globular head; SI) from a consensus among PEDV non S-INDEL strains, 5) SI from a 

543 consensus among PEDV S-INDEL strains. 

544 Figure 2. PEDV multiplex (6-plex) fluorescent microbead-based immunoassay (FMIA) sample- 

545 to-positive (S/P) ratios of serum antibody (IgG) response (mean, SE) over time against the 5 

546 recombinant polypeptides (SI non S-INDEL, SI S-INDEL, N, M, and E), and the whole virus 

547 (WV) antigen in pigs (n = 12 per group) inoculated with PEDV (USA/IN/2013/19338E), TGEV 

548 Miller (ATCC VR-1740), TGEV Purdue (ATCC VR-763), PRCV (ATCC VR-2384), PDCoV 

549 (USA/IL/2014), and a negative control (sham inoculated). Different letters denoted statistical 

550 differences (p=< 0.05). Samples above the estimated S/P cut-off (dashed line) were considered 

551 positive. 

552 Figure 3. Distribution of cumulative FMIA IgG sample-to-positive (S/P) ratios obtained for each 

553 recombinant polypeptide (SI non S-INDEL, SI S-INDEL, N, M, E) on serum samples (n = 792 

554 total; n = 132 per group) collected at DPI -7, 0, 3, 7, 10, 14, 17, 21, 28, 35, and 42 from pigs (n = 

555 72 total; n = 12 per group) inoculated with PEDV (USA/IN/2013/19338E), TGEV Miller (ATCC 

556 VR-1740), TGEV Purdue (ATCC VR-763), PRCV (ATCC VR-2384), PDCoV (USA/IL/2014), 

557 and a negative control (sham inoculated). The FMIA S/P cut-off values estimated for each 

558 individual antigen are presented in the graph (dashed line). 
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559 Figure 4. Performance of PEDV whole virus (WV) indirect ELISA on experimental samples of 

560 precisely known porcine coronavirus infectious status, i.e., pigs inoculated with PEDV 

561 (USA/IN/2013/19338E; n = 12), TGEV Miller (ATCC VR-1740; n = 12), TGEV Purdue (ATCC 

562 VR-763; n = 12), PRCV (ATCC VR-2384; n = 12), PDCoV (USA/IL/2014; n = 12), and a 

563 negative control (sham inoculated; n = 12). A) Sample-to-positive (S/P) ratios of serum antibody 

564 (IgG) response (mean, SE) over time in each inoculation group. B) Distribution of cumulative 

565 ELISA WV IgG sample-to-positive (S/P) on serum samples (n = 792 total; n = 132 per group) 

566 collected at DPI -7, 0, 3, 7, 10, 14, 17, 21, 28, 35, and 42. Different letters denoted statistical 

567 differences (p=^ 0.05). Samples above the estimated S/P cut-off (dashed line) were considered 

568 positive. 

569 




Downloaded from http://jcm.asm.org/on February 17, 2017 by UNIV OF CALIF SAN DIEGO 



Figure 1 


1. PEDVrSI S-INDEL 

2. PEDV rSI non S-INDEL 

3. PEDV rN 

4. PEDV rM 

5. PEDV rE 


1. PEDVrSI S-INDEL 

2. PEDV rSI non S-INDEL 

3. PEDV rN 

4. PEDV rM 


1. PEDVrSI S-INDEL 

2. PEDV rSI non S-INDEL 

3. PEDV rN 


1. PEDV rSI S-INDEL 

2. PEDV rSI non S-INDEL 

3. PEDV rN 

1. PEDV rSI S-INDEL 

2. PEDV rSI non S-INDEL 

3. PEDV rN 


1. PEDVrSI S-INDEL 

2. PEDV rSI non S-INDEL 


1. PEDV rSI S-INDEL 

2. PEDV rSI non S-INDEL 


1. PEDVrSI S-INDEL 

2. PEDV rSI non S-INDEL 


1 10 20 30 40 SO 60 70 80 90 

QDVTRCQSTi NFRRFFSKFN VQAPAVVVLG GYLPSMNSSS WVCGTGL ETA SGVHGIFLSY IDAGQGFEIG ISQEPFDPSG YQLYLHKATN 

QDVTRC SANT NFRRFFSKFN VQAPAVVVLG GYLPIGENQG VNSTWYCAGQ HPTASGVHGI FVSHIRGGHG FEIGISQEPF DPSGVQLYLH 

MASVS FQDRG RKRVPLSLYA PLRVTNDKPL SKVLANNAVP TNKGNKDQQI GYWNEQIRWR MRRGER I EQP SNWHFYYLGT GPHADLRYRT 

FVNSIRLWRR THSWWSFNPE TDALLTTSVM GRQVCIPVLG APTGVTLTL L SGTLLVEGYK VATGVQVSQL PNFVTVAKAT TTIVYGRVGR 

MLQLVNDNGL VVNVILWLFV LFFLLIISIT FVQLVNLCFT CHRLCNSAVY TPIGRLYRVY KSYMRIDPLP STVIDV 

100 110 120 130 140 ISO 160 170 1B0 

GNHNAIARLR ICQFPDNKTL G PTVNDVTT G RNCLFNKAIP AYMQDGKNI V VGITWDNDRV TV F AD K IYH F YLKNDWSRVA TRCYNKR SCA 

K ATNGNTNAT ARLR ICQFPS IKTLGPTANN DVTTGRNCLF NKAIPAHMSE HSVVGITWDN DRVTVF SDK I YYFYFKNDWS RVATKCYNSG 

RTEGVFWVAK EGAKTEPTNL GVRKASEKP I IPNFSQQLPS VVEIVEPNTP PTSRANSRSR SRGNGNNRSR SPSNNRGNNQ S RGNSQNRGN 

SVNASSGTGW AFYVRSKHGD YSAVSNPSSV LTDSEKVLHL V 

190 200 210 220 230 240 250 260 270 

MQYVYTPTYY MLNVTSAGED GIYYEPCTAN C S GYAANVF A TD S NGHI PEG FSFNNWFLLS NDSTLLHGKV VSNQPLLVNC LLAIPKIYGL 

GCAMQYVYEP TYYML NVTSA GEDGISYQPC TANCIGYAAN VFATEPNGHI PEGFSFNNWF LLSNDSTLVH GKVVSNQP L L VNCLLAIPKI 

NQGRGASQNR GGNNNNNNK S RNQSKNRNQS NDRGGVTS RD DLVAAVKDAL KSLGIGENPD KLKQQQKPKQ ERSDSSGK NT PKKNKSRATS 

280 290 300 310 320 330 340 350 360 

GQFFSFNQTM DGVCNGAAAQ RAP EAL R F Ni NDTSVILAEG S IVLHTALGT NLSFVCSNSS DPHLATFAIP LGATQVPYYC FLKVDTYNST 

YGLGQFFSFN QTI DGVCNGA AVQRAPEALR FNINDTSVIL AEGS IVLHTA LGTNFSFVCS NSSNPHLATF A I P L GATQV P YYCFFKVDTY 

KERDLKDI PE WRRIPKGENS VAACFGPRGG FKNFGDAEFV E KGVDASGYA QI AS LAP NVA ALLFGGNVAV RELADSYEIT YNYKMTVP K S 

370 380 390 400 410 420 430 440 450 

VYKFLAVLPP TVREIVITKY GDVYVNGFGY LHLGLLDAVT I N FTGHGTDD DVSGFWTIAS TNFVDALIEV QGTAIQR I LY CDDPVSQLKC 

N STVYK F LAV LPPTVREIVI TK YGDVYVNG FGYLHLGLLD AVTINFTGHG TDDDVSGFWT IASTNFVDAL IEVQGTAIQR ILYCDDPVSQ 

DPNVELLVSQ VDAF KTGNAK PQRKKEKKNK RETTQQLNEE AIYDDVGVPS DVTHANLEWD T AVDGGDT AV EIINEIFDTG N 

460 470 480 490 500 510 520 530 540 

SQVAFDLDDG FYP I SSRNLL SHEQPISFVT LPSFNDHSFV NI TVSAS F GG H S GAN L I AS 6 TTINGFSSFC VDTRQFTISL F YNVTNSYGY 

LKCSQVAFDL DDGFYPISSR NLLSHEQPIS FVTLPSFNDH SFVNITVSAS FGGHSGANLI ASDTTINGFS SFCVDTRQFT ISLFYNVTNS 

550 560 570 580 590 600 610 620 630 

VSKSQDSNCP FTLQSVNDYL SFSKFCVSTS L LASACT I D L FGYPEFGSGV KFTSLYFQFT KGEL ITGTPK PLEGVTDVSF MTLDVCTKYT 

YGYVSKSQDS NCPFTLQSVN DYLSFSKFCV STS L LASACT I DLFGYPEFG SGVKFTSLYF QFTKGELITG TPKPLEGVTD VSFMTLDVCT 

640 650 660 670 680 690 700 710 717 

lYGFKGEGli TLTNSSFLAG VYYTSDSGQL LAFKNVTSGD VYSVTPCSFS E QAAYVDDD \ VGVISSLSSS TFNSTRELPG FFYH 

KYTIYGFKGE GIITLTNSSF LAGVYYTSDS GQLLAFKNVT SGAVYSVTPC SFSEQAAYVD DDIVGVISSL SSSTFNSTRE LPGFFYH 


Downloaded from http://jcm.asm.org/on February 17, 2017 by UNIV OF CALIF SAN DIEGO 



Figure 2 


PEDV SI (S-INDEL) 



PEDV SI (non S-INDEL) 




PEDV 

TGEV-Purdue 
TGEV-Miller 
PRCV 
PDCoV 

Negative control 


Day post Inoculation 


Day post Inoculation 


Day post inoculation 


PEDV WV particles 



Day post Inoculation 


H"*- i + i l i i H 


3 7 10 14 17 21 28 35 42 

Day post inoculation 



3 7 10 14 17 21 28 35 42 

Day post Inoculation 


Downloaded from http://jcm.asm.org/on February 17, 2017 by UNIV OF CALIF SAN DIEGO 






































Figure 3 
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Figure 4 
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